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We analyze experimentally chemical waves propagation in the disordered flow field of a porous 
medium. The reaction fronts travel at a constant velocity which drastically depends on the mean flow 
direction and rate. The fronts may propagate either downstream and upstream but, surprisingly, 
they remain static over a range of flow rate values. Resulting from the competition between the 
chemical reaction and the disordered flow field, these frozen fronts display a particular sawtooth 
shape. The frozen regime is likely to be associated with front pinning in low velocity zones, the 
number of which varies with the ratio of the mean flow and the chemical front velocities. 
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Introduction. Interface motion and front propagation 
are relevant to many processes, including chemical reac- 
tions Q , population dynamics in biology, flame propaga- 
tion in combustion Q and chemotaxis [Bj]. The dynamic 
of front propagation is well understood in stagnant fluids 
@, Q , whereas the effect of fluid flow on front motion is 
currently of immense interest 0-0. In this experiment, 
we consider an autocatalytic reaction between two chemi- 
cal species which generates a self-sustained reaction front 
[13, LLL|. These fronts propagate as solitary waves with 
a constant velocity (V x oc ^/D m a) and a stationary con- 
centration profile of width (l x = D m /V x ), reflecting both 
the balance between molecular diffusion D m and reaction 
rate a. 

In heterogenous flow field, non-linear coupling between 
reaction, diffusion and advection generally induce long- 
range deformation of the reaction front. For instance, in 
a parabolic Poiseuille flow 0, EH, [l3j , this coupling leads 
to a stationary front concentration profile propagating at 
a constant velocity, which shape and velocity are related 
to the local flow structure and the boundary conditions. 
More complex flows have been addressed in porous media 
[Hi [H| or in cellular flow \^ . Interestingly, both these 
experiments have shown the existence of static fronts for 
flow in direction opposite to the chemical wave propaga- 
tion. In this context, we have investigated the interac- 
tion of a self-sustained reaction front with the disordered 
flow field in a porous medium. Our experiments reveal 
new features of front propagation in porous media due to 
the use of a wide and transparent quasi-two dimensional 
design, allowing the flow velocity fluctuations measure- 
ments and the direct observation of the reaction fronts. 
Depending on the flow velocity and direction, the self- 
sustained fronts can travel downstream, upstream with 
a constant front velocity but they can also remain static 
over a range of flow rate values. 

In this Letter, we show visual aspects of frozen re- 
action fronts in porous media. These static fronts are 
strongly distorted and display a remarkably straight saw- 
tooth pattern, as shown in Fig. [TJ the patterns aspect, 
the number of peaks and size, evolving with the mean 



flow rate. We determine first the main characteristics of 
the flow field such as its correlation length and the spatio- 
temporal maps of the local velocities along the front as 
well as its probability distribution function (PDF). Ana- 
lyzing the freezing dynamics then provides a better phys- 
ical understanding of frozen fronts formation mechanism 
in porous media. 
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Figure 1: Typical self-sustained chemical fronts in porous me- 
dia. The first label corresponds to the mean flow direction, 
supportive S or adverse A to the chemical wave propagation; 
the second label denotes the front propagation, downstream 
D, upstream U or static S. 



Experimental setup. We performed experiments with 
the lodate Arsenous Acid (IAA) autocatalytic reaction 
in iodide (I"): 3H 3 As0 3 + IO3 + 51" — > 3H 3 As0 4 + 
6I~. For the concentration values used here ([IO3 ]o = 7.5 
mM, [HsAsOaJo = 25 mM), the arsenous is in excess 
The front position is localized by the transient generation 
of iodine during the reaction and marked by polyvinyl 
alcohol (see Fig [T] for instance) [l6| . 

The porous medium consists of 50% mixture of 1.5 
and 2 mm glass spheres (porosity: 48 ± 2%). They are 
packed in a transparent rectangular cell (30 x 10 x 1 cm 3 ) 
as shown in Fig. |2J To overcome buoyancy effect [I?}, we 
increase the viscosity to 10 mPa.s. The cell is filled with 
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Figure 2: (a) Experimental set-up schematic, (b) Spatio- 
temporal map of the velocity variations measured with dye 
transport in the porous medium, (c) PDF of the above ve- 
locities map (dots) and log- normal fit (dashed line). 



reactant solution and the chemical reaction is initiated 
at its bottom end. In the absence of flow, the reaction 
develops into a flat horizontal front propagating upward 
at the chemical velocity V x ~ 10 fim/s and l x ~ 100 fim. 
When the front reaches the desired vertical position, the 
flow is switched on and the fluid can be either sucked 
out or injected parallel to the vertical from the top of the 
cell. We extract the location h(x) of the chemical pseudo- 
interface at a given transverse distance x and define the 
front velocity by v(x,t) = (h(x,t + At) — h(x,t))/ At; 
At is such that the variation of h(x) is larger than the 
correlation length. 

The net flow is characterized by its average value, U 
and the characteristic length of the velocity fluctuation 
Id, is inferred from tracer dispersion experiments [l8[[l9|: 
Id = (1.8 ± 0.1) mm which is close to the mean particle 
diameter, in agreement with [20j |. Tracer front velocity 
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Figure 3: Front velocity Vf versus mean flow velocity in the 
porous medium, u. 

fluctuations u(x,t) are shown in Fig. [2j The spatiotem- 
poral map displays the flow heterogeneities of character- 
istic size larger than Id and the flow velocities PDF is 
asymmetric and well fitted by a log-normal distribution 
(solid line). 

Front velocity. We performed experiments over a wide 



range of mean flow velocities and measured the reac- 
tion front velocity. The fronts propagate with a veloc- 
ity slightly fluctuating around a mean value, noted Vf. 
Counting as positive velocities in the direction of the re- 
action front propagation, U > corresponds to support- 
ive flow (S) and U < to adverse flow (A). The mean 
front and mean flow velocities are, now, normalized by 
V x ; u = U/V x and Vf = VfjV x . As can be seen from 
Fig. [3j Vf depends on both the value and direction of u 
and one can identify four different regimes of front prop- 
agation; the corresponding front's aspects are shown in 
Fig. [TJ For supportive flow (u > 0), the front propa- 
gates downstream, i.e. in the flow direction (vf > 0). 
In this supportive downstream regime (SD) , Vf varies 
almost linearly with u with a slope larger than one. As 
demonstrated in Poiseuille flow 0, [l2| , this may be at- 
tributed to the selection of the fastest streamline by the 
reaction front when propagating downstream. For ad- 
verse flow (u < 0), the competition between the chemical 
wave and the flow field leads to more complex behaviors. 
For a small adverse flow rate — 1 < u < 0, the chemical 
wave overcomes the flow leading to upstream propagat- 
ing fronts (AU Vf > 0); Vf varies smoothly across the 
regimes SD and AU and approaches zero. The most 
salient regime is observed below u c± —1. The reac- 
tion fronts are static without any fluctuations and dis- 
play frozen sawtooth like patterns (Fig. [1]). They are ob- 
served over a range of mean flow velocities u G [—3.2,-1] 
which corresponds to the plateau domain (vf = 0, AS 
regime). Finally, for u < —3.2, the flow is strong enough 
to overcome the chemical wave and so that the front 
moves downstream, opposite to the chemical wave veloc- 
ity (vf < 0, AD). Its velocity varies again almost linearly 
with u. However, in contrast to the SD regime, the re- 
sulting slope is smaller than one, which is in agreement 
with similar experiments using Hele-Shaw cells [7|. 
Static fronts have been observed in Poiseuille flow in 
Hele-Shaw cells @, G3, but they were obtained only for 
a unique value of u. In contrast, the velocity diagram 
determined here is similar to that observed in other con- 
texts. For instance, the domains of observation of frozen 
fronts has been determined in combined cellular and 
mean opposite flow [9[ as a function of the vortex and 
the mean adverse flow intensities. Static fronts in packed- 
bed reactors with different chemistries [3, [H| have also 
been reported for adverse flow. Quite remarkably the 
reaction fronts exhibit similar behaviors despite the im- 
portant differences between these systems. Furthermore, 
one common feature of all these experiments is the het- 
erogeneity of the flow and shows here the key part the 
flow spatial structures play in frozen fronts formation. 

Spatio-temporal dynamics. One of the major motiva- 
tions of these experiments, is to understand how the flow 
velocity fluctuations alter the front propagation and ac- 
count for the different regimes. This can be inferred from 
the fronts dynamics; Fig. H] displays the spatio-temporal 
fluctuations of v(x,i) and the associated PDFs in the 
different regimes. For SD fronts, the mean flow and the 
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Figure 4: Spatiotemporal representation of v(x,t) for the four regimes. Color code intensity represents the magnitude of the 
velocity; red: v > 0, blue: v < 0. White regions corresponds to v — 0. The corresponding PDFs of v(x,t) (dots) and log-normal 
fit (dashed line). 



chemical wave are in the same direction, therefore the 
measured local front velocities, v(x,t), have positive val- 
ues during the experiment. The PDF of these velocities 
has a log-normal shape with almost the same standard 
deviation as the one measured previously for the hydro- 
dynamic flow field. This suggests that the hydrodynamic 
fluctuations dominate the front dynamic in this regime. 
The reaction front propagation is strikingly different in 
adverse flow regimes. The PDFs are modified indicating 
a more complex dynamical response than just the addi- 
tion of the chemical velocity to the flow velocity distri- 
bution. The spatio-temporal maps exhibit zero velocity 
regions. This is a signature of transient static portions of 
the fronts, resulting in a local maximum at v = on the 
PDFs. In the AU regime, one can identify long range 
lateral correlation of the velocities as dark regions on 
the spatio-temporal. They indicate transversal propaga- 
tion of the reaction at constant rate. At a given location, 
when the front stops, it remains static until laterally trav- 
eling parts of the front reach these waiting regions. In 
this configuration, the fronts propagate upstream, how- 
ever the PDF displays a few negative velocities. Those 
values can be attributed to the locations where the flow 
velocity is locally larger than the chemical wave veloc- 
ity. In the AS regime, one can note that the final static 
pattern is achieved after a transient phase. Therefore, 
the spatio-temporal map describes the freezing dynam- 
ics of the front. Velocity fluctuations are negative, il- 
lustrating the fact that the front is mostly pushed back 
by the flow until it pins locally to a point-like region (at 
X — 63 mm in Fig. HJ. The final front structure displays 
a sawtooth shape (Fig. [1]), where each pinning point is 
at the vertex of a peak. However, before achieving the 
final pattern, transitorily pinned parts of the front are 
visible on the map as white stripes. In contrast to the 
AU regime, this demonstrates the long transversal corre- 
lation of zero velocity regions over time. Hence, the front 



develops transient static parts during the whole experi- 
ment until it reaches an entirely frozen state. Finally, 
in the AD regime, fronts move downstream all along the 
experiments, and thus in this regime the flow seems to 
dominate the propagation of reaction. Yet a small peak 
at v = on the PDF indicates the persistence of tran- 
sitorily pinned parts of the fronts, which progressively 
vanish as the mean adverse flow velocity is increased. 
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Figure 5: Left: Frozen fronts aspects for different mean flow 
rate: from top to bottom, u = —1.4, —1.6, —1.9, —2.5, —3.2. 
Right: Measurement of versus u; solid line, theory. 



Discussion. The dynamical study has revealed pin- 
ning regions of the reaction fronts in adverse flows. 
They play a determining role in the freezing mecha- 
nism for the plateau regime (AS). Experiments on sta- 
tionary Belousov-Zhabotinsky pulses have been reported 
in packed bead reactors Q an d linked to an "excited 
stagnant pockets" mechanism. When the reaction front 
reaches these regions, they act as point source of chem- 
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ical reaction. Producing and spreading catalyst, these 
pockets can sustain the reaction front in their vicinity. 
Although the chemical reaction is very different in that 
case, one can follow a similar argument in our experi- 
ments such that these pinning zones may correspond to 
stagnant flow pockets. 

However this argument cannot account for reaction 
fronts moving downstream for sufficiently large adverse 
flow. This indicates that pinning points must be asso- 
ciated with weak but generally non-zero velocity zones, 
i.e. —U(f) < V x . Depinning of the front then occurs 
when the local flow velocity reaches a critical value which 
pushes the reaction out of this region. In this case, in- 
creasing the mean velocity results in reducing the number 
of low velocity zones required for pinning. This mecha- 
nism is supported by visual observations displayed on 
Fig. [5] Several frozen patterns recorded for different 
u values in the plateau domain are plotted and one can 
identify the location of the pinning points with the vertex 
of the front's peaks. Note that the two lateral boundaries 
can also be considered as pinning zones. As the flow ve- 
locity increases, the number of pinning points decreases 
until they vanish completely (except the two cell lateral 
boundaries) for u = —3.2. This suggests that, statisti- 
cally, the chemical reaction finds fewer pockets which can 
remain excited when u increases. The shapes of the final 
static fronts depend on the distribution of these pinning 
points but also on the initial location of the front. 

Away from these zones, the flow velocity is much higher 
than V x , and yet, the fronts reach a frozen state all over 
the porous medium. In order to understand these frozen 
fronts, one has to deal with the two flow regimes already 
analyzed in simple flows 0, G3 . In the so called mixing 
regime, the size of the flow extension (Id) is much smaller 
than the chemical front width, Z x , molecular diffusion has 
enough time to transversally mix the reactants leading to 
the Galilean invar iance: Vf = U + V x . Therefore, static 
fronts are possible for a unique flow rate U = — V x ; even 
with the flow distribution of Fig. [2j it is not enough to 
account for the plateau. In the opposite, (Id » l x ), thin 
front eikonal limit [21], at each point of the front surface, 
the normal component of the local interface velocity sat- 
isfies : 

V f (r) -n = V x + U(r)-n + D m K (1) 

where U(r) is the local fluid velocity at the front posi- 
tion r, n is the local unit vector normal to the interface 
and ft, the curvature of the interface. In this regime, 
frozen fronts can adjust their local inclination and cur- 



vature to accommodate larger local flow velocities and 
satisfy: Vf = in Eq. [TJ Hence, increasing the flow 
intensity from u = —1 to u = —3.2 results in larger ad- 
verse flow leading to more inclination and bending of the 
front around pinning zones. Note that, we are here in an 
intermediate regime (Id ~ Wl x ) where reactants concen- 
trations variations occur over a few / x , however the thin 
front limit still explains the aspect of the frozen states. 
In Fig. [5j the teeth angle decreases with \u\. Its value 
is in agreement with: V x + U sin(0/2) = (solid line), 
which derives from Eq. [T] with the curvature term ne- 
glected on the straight parts of the front. Particularly, 
this pattern is analogous to the V-shape displayed on 
flame fronts, strongly inclined in order to counterbalance 
high flow rate 

Finally, we have shown that the transient static por- 
tions of the fronts observed on the spatio-temporal maps, 
correspond to transitory pinning zones. During the ex- 
periments, the reaction fronts are temporally pinned at 
these locations. They form local static parts without 
leading to a completely frozen pattern and depin sud- 
denly causing burst-like events. The irregular propaga- 
tion of the fronts in AS and AD regimes is hence due to 
successive rearrangements of the fronts when they travel 
from one low velocity zone to another one. In these spe- 
cific regimes, the fronts propagate through local pinning 
and depinning events whose statistical properties [23j de- 
serve further investigations aided by numerical simula- 
tions. 

Conclusion. We have experimentally investigated the 
propagation of chemical waves coupled to the disordered 
flow of a model porous medium. Depending on the mean 
flow direction and intensity, the reaction fronts exhibit 
four different regimes. Our experiments show the key 
role played by the flow heterogeneities on the chemical 
front dynamics. An interesting freezing regime is ob- 
served and explained by front pinning to low flow velocity 
regions. More specifically, we have analyzed the complex 
dynamic features involved in this pinning process. We 
thus expect that the domain of existence of this regime 
strongly depends on the flow heterogeneities distribution 
and further numerical simulations should be performed 
to analyze this dependence. 
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